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1 Introduction 



The analysis of singular systems is an interesting problem in Physics, as such systems 
appear in many relevant physical problems. Such analysis is usually carried out using the 
generalized Hamiltonian formulation, developed by Dirac [1, 2, 3, 4], where the canonical 
Hamiltonian is not uniquely determined due to the singularity of the Hessian matrix; what 
causes the appearance of relations between canonical variables. These constraints, mul- 
tiplied by Lagrange multipliers, are added to the canonical Hamiltonian and consistency 
conditions are implemented to eliminate some degrees of freedom of the system. 

Despite the outstanding success of Dirac's formalism, new approaches to the analysis 
of singular systems are always welcome because they may reveal new mathematical and 
physical information about the system in study. Among others, an alternative method 
to analyze singular systems is the Hamilton- Jacob i formalism [5, 6], which has been used 
in many examples [7, 8, 9] and generalized to higher order singular systems [10, 11] and 
systems with Berezinian variables [12]. This formalism uses Caratheodory's equivalent 
Lagrangians method [13] to write a set of Hamilton- Jacobi partial differential equations 
from which one can obtain the equations of motion as a set of total differential equations 
in many variables. 

One example of physical system described by a singular Lagrangian, and that has 
already been studied through Dirac's method [14, 15, 16], is the Teleparallel Equivalent 
of General Relativity (TEGR) which is an alternative formulation of General Relativity 
[17] developed in Weitzenbock space-time [18]. In oposition to General Relativity, in 
TEGR the curvature tensor vanishes but the torsion tensor does not so that, in this 
geometrical framework, the gravitational effects are caused by the torsion tensor and not 
by curvature. 
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The TEGR has been successfully analyzed through Dirac's Hamiltonian formalism 
[14, 15] and generated successfull applications [19]. Moreover, many of the characteristics 
regarding the interaction of spin 0, 1 and spinor fields in TEGR have been studied recently 
[20, 21, 22, 23, 24] as well as its gauge symmetries [25]. Others aspects of TEGR, as the 
energy-momentum tensor and geodesies and "force" equation are addressed in reference 
[26] and in the references cited therein. 

Our intention in this work is to add a different point of view to the analysis of TEGR by 
studying it through the above mentioned Hamilton- Jacobi formalism for singular sistems. 
First, we will introduce the Lagrangian density of TEGR in a form which is appropri- 
ate to our approach. Then we address the basic aspects of Hamilton- Jacobi formalism 
for singular systems and apply such formalism to TEGR. Finally, we present our final 
comments. 



2 The Lagrangian of TEGR 

In this section we summarize the Lagrangian formulation of TEGR in terms of the tetrad 
field, as presented in reference [15] where a global 5*0(3, 1) symmetry is taken from outset. 
This choice is done because, when starting from a local SO (3, 1) symmetry, it may not be 
possible, with certainty reference systems choices, to obtain a set of first class constraints 
[14]. 

So, we take the Lagrangian density of TEGR in empty space-time given, in terms of 
the tetrad field e afi , by 

L(e) = -keJ: abc T abc , (1) 

where Latin letters are SO(3, 1) indexes (taking values (0), ...(3)), Greek letters are space- 
time ones (taking values 0, ...,3), e = det(e a M ), k = and G is the gravitational 
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constant. Besides that, the torsion tensor T abc = eb^e c u T afMU is defined in terms of the 
tetrad field as 

and its trace is defined as 

Ti rpa 
b — J- ab, 

while the tensor S abc is defined as 

Y^abc \_ ^rpabc _|_ rpbac rpcab^ _|_ ^_ ^q acr J^ } Tj a ^T C ^J 
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such that 



Sabcrp 1 rpabcrp , rpabcrp rparp 

J-abc — l-abc-r^ 1 bac~ 1 1 a- 



The fields equations can be obtained from the variation of L with respect to e afl and 
are equivalent to Einstein's equations in tetrad form [14] 
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= lej^^-Ie^e)}. (2) 



Let us now consider the tetrad field e a/1 in terms of the 3 + 1 decomposition. In this 
case the space-time manifold is assumed to be topologically equivalent to M x R, where 
M is a noncompact three-dimensional manifold. We consider that 4 e a ^ is a tetrad field 
for M x R. In terms of the lapso iV and shift N l (i, j, k... = 1, 3) functions we have 



V , = e° 
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4 e ai = e ai + (jV*/JV)77°, 

V = e a i N i + r ] a N, 



e m = e\g ki , 



rf = -N 4 e a0 , 



where e ai and e m are now restricted to M. Moreover, rf is an unit timelike vector such 
that r] a e a i = and g lk is the inverse of g ik = e ai e a k . The determinant of the tetrad field 
is now given by 4 e = A^e, where e = det(e® k ). In terms of the 3+1 decomposition, the 
TEGR Lagrangian density L can be written as 



L = W + e ( ° )l e (0)fc / (0 il {l)k + e^e^ k l (n)l l m - 2e^% yi l (n)k ) 

-AkeY: a ^H ai - kNeZ abc T abc 
where we made the use of the following definitions 

Y^abc ^_ ^rpabc _|_ rpbac rpcab^ _|_ ^_ ^p acr J } b Tj'^T'^j 

rpabc _ g^^g c ^^~' ct 

The momenta Tl ak conjugated to e ak are obtained by 

U ak = 5L 

8e ak 

where 
end 

n^ s = ~ (l^ s + e (0)l e (0)s / (r) i + e {r)% e {n)s l {n)t - 2/e (r)s ) + 4A;e£ (r)(0)s . 

Note that here fi* = g ik l^ k and I = e^H {l)i . 

In order to simplify the calculations we can make a choice of reference frame, analogous 
to the one made in ADM formulation. This choice is usually refered in literature as the 
Schwinger's time gauge [28], e^i = e( ) 1 = 0, so that the momenta defined above are 



such that U ks = Tl sk and n(°) fc = ^ k . In this case it is possible to write the Lagrangian 
density L as [14] 

L = U^e {l)i - H c , (3) 

where 

H c = NC + N'Q, (4) 
C = ± (u^IUj ~ ^n 2 ) + eY? ik T ijk - 2d t (eT), (5) 

and 

d = e {l)i d k U^ k + U^ k T m . (6) 

In Eq. (3) there is no time derivatives of the functions iV and N l , so the Lagrangian 
density L is singular. With these results, we now can investigate the integrability condi- 
tions of the Lagrangian density of the TERG, given in Eq. (3), in the Hamilton- Jacobi 
formalism, what it will be done in the next section. 



3 The Hamilton-Jacobi Formalism 

The Hamilton-Jacobi formalism, recently developed to analyze singular systems [5, 6, 
10, 12], uses the equivalent Lagrangians method [13] to obtain a set of Hamilton-Jacobi 
partial differential equations [5, 6]. We suggest the references just mentioned for details 
and present here only the main aspects of this formalism. For this, let us consider a 
singular Lagrangian function L = L(qi,<ji,t), where % = 1,...,N. The Hessian matrix is 
then given by 

Being the rank of the Hessian matrix P = N — R < N, we can define, without loss of 
generalization, the order of the variables qi in a such way that the P x P matrix in the 



right bottom corner of the Hessian matrix be nonsingular. So there will be R relations 
among canonical variables given by 

p a = -H a (V; p a ) ; a = 1, ...,R; (8) 

which correspond to Dirac's primary constraints $ a = p a + H a (q l ; p a ) ps 0. From these 
we get the Hamilton- Jacobi partial differential equations, given by 

H' a = H a (t, qi ,p a )+p a = 0, (9) 

H' = H c {t,q u p a )+p = 0, (fO) 

where 

dS(t,qj) 

Po = (12) 

and i, j = 1, N, a — 1, i?; a = i? + 1, iV; if c is the canonical Hamiltonian and S 
is the Hamilton principal function. 

It can be shown that the equations of motion are total differential equations for the 
characteristics curves of the differential partial equations (9) and (10), being given by [5] 

dq% = ^f dt + W dq °" (13) 

a tti a tti 

d Pt = -^dt-^ dqa , (14) 

where for i = 1,...,R equation (13) above becomes a trivial identity. Using standard 
techniques of partial differential equations, it can be shown [12] that the equations above 
are integrable if and only if the functions H' a satisfy 

dH' p = {H' p ,H' a }dt a , (15) 
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where a, (5 — 0, 1, ...,R; t — t (so that t v = (t,q a )) and the symbol {..., ...} denotes the 
Poisson bracket defined on the phase space of 2N + 2 dimension that includes t — t and 
its canonical momentum po. 



3.1 The Teleparallel Lagrangian case 



Let us now consider the Lagrangian density given in Eq. (3) in the context of the 
Hamilton- Jacobi formalism. For this purpose we note that this Lagrangian density does 
not depend on time derivatives of the shift and lapso functions, therefore we define the 
set of Hamilton- Jacobi partial differential equations as 



H' 0c = I d 3 x(H c (x)+ Po (x)) = 0, 



(16) 



H' Q = J d 3 xTl (x) = 0, (17) 

H'i = J d 3 xUi(x) = 0, i = l, ..3; (18) 

where n and Hi are the canonical momenta conjugated to the shift and lapso functions, 
respectively; H c is the canonical Hamiltonian andpo is the "momentum" conjugated to the 
time parameter. So, in this approach, both iV and iV* are taken as evolution parameters, 
together with t, from the beginning. Now we apply the integrability conditions given by 
Eq. (15) to the equations above. After some calculation we obtain 

(19) 
(20) 

dt = 0, (21) 
(22) 



dH' = {Hq, H' 0c } dt = -Q d 3 xC(x)j dt = 0, 



dH[ = {Hi, H' 0c } dt = -(J d 3 xCi(x^j dt = 0, 
dH' 0c = {H' 0c , H' 0c } dt + {H' c , H' } dN + {H' c , H>} dN* = 
where 



d : '\r— ( : r 
ox 



I = N 2 { -5% n - L ,r „„ ■ T„„ ' ) II'"'. m.j I. ^ 3: 



with the brackets on the indexes indicating antisymmetrization. Equation (21) above can 
be transformed in a surface integral and, therefore, at the surface of integration we must 
have 



However, the latest equation is a consequence of the fact that the canonical momenta 
jjmj j n {.jjis approach are symmetrical, so the antisymmetric components n' m ^ must van- 
ish, therefore the Equation (23) is satisfied in in the whole spacetime, what corresponds 
to primary constraints in the approach of reference [14]. 

The integrability conditions given by equations (19), (20) and (21) imply that C — 0, 
d = and n^'l = and no new conditions arise. The integrability conditions are 
equivalent to the consistency conditions obtained in reference [14] using Dirac's method, 
and the quantities C, Cj, U\ m ^ constitute a set of first class constraints. 

4 Final remarks 

In this work we analyzed the Lagrangian density of TEGR, with a specific choice of ref- 
erence frame, by using the Hamilton- Jacobi formalism, which was recently developed to 
treat singular systems. Such Lagrangian has already been studied through Dirac's Hamil- 
tonian formalism, where the consistency conditions produce a set of first class constraints 
[14]. In our analysis, the integrability conditions of Hamilton- Jacobi formalism produce 
results that are identical to those obtained in reference [14] through the use Dirac's Hamil- 
tonian method. 

However, one of the most interesting characteristics of Hamilton- Jacobi formalism is 
the possibility of avoiding specific choices of gauge or reference systems. So, we believe 
to be possible to study the Lagrangian density of TEGR without assuming any a priori 




(23) 
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restriction on the tetrad fields, like Schwinger's time gauge. Our expectation is that such 
restrictions should naturally arise as consequence of integrability conditions in Hamilton- 
Jacobi formalism, as happens in other singular systems [29, 30]. This question is presently 
under our study. 
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